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Abstract 

We examine the charm component F2^c{x., Q^, m^) of the proton structure function 
F2(x,Q^) in three different schemes and compare the results with the data in the x 
and region explored by the HERA experiments. Studied are (1) the three flavour 
number scheme (TENS) where the production mechanisms are given by the photon- 
gluon fusion process and the higher order reactions with three light-flavour parton 
densities as input (2) the four flavour number scheme (EENS) where F2^c is expressed 
in four light flavour densities including one for the charm quark and (3) a variable- 
flavour number scheme (VENS) which interpolates between the latter two. Both the 
VENS and the TENS give good descriptions of the experimental data. However one 
cannot use the EENS for the description of the data at small Q^. 
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The study of charm production in deep- inelastic electron-proton scattering has 
become an important issue in the extraction of parton densities in the proton. The 
reason is that the charm content F2^c{x,Q'^ j^"^) {m is the mass of the charm quark) 
in the proton structure function ^2(2;, Q^) has grown from around one percent in the 
X and region covered by the EMC experiment [|T| to around twenty-five percent 
in the x and region covered by the HI and ZEUS [Q experiments at HERA. 
Therefore the analysis of data to yield parton densities can no longer treat charm 
electroproduction as a small correction. 

Let us begin with a brief review of the QCD contributions to F2^c{x,Q'^,m^)- In 
lowest order (LO) the charm quark pair is produced via the photon-gluon fusion 
mechanism which implies that the Born contribution only contains a gluon in the ini- 
tial state (extrinsic production). In next-to- leading order (NLO) we have in addition 
to the above reaction also Bethe-Heitler and Compton subprocesses containing the 
three light flavours u, d, s and their antiparticles u, d, s in the initial state. We will 
call this the three flavour number scheme (TENS) description for charm production. 
The NLO calculations |^ reveal that the photon-gluon process dominates the other 
production mechanisms so that charm electroproduction yields a measurement of the 
gluon density, in particular at small x. Generalizing the TENS to all orders in we 
obtain the following formula for the charm content of the proton structure function 
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where stands for the charge of the charm quark and the number of light flavours 
is three (nj = 3). The summation variable k refers to the light quarks u,d and s. 
The variable z is the partonic longitudinal momentum fraction which has a maximum 
value Zmax = Q"^ /{Q"^ + 4m^). The function A is the non-singlet (NS) (with respect 
to the flavour group) combination of light flavour densities, which are functions of 
the scale /x. The function S is the singlet (S) combination of these densities while 
G stands for the density of the gluon. Eurther L2,fc and H2^k{k = q,g) represent 
the charm-quark coefficient functions which can also be separated into flavour singlet 
and flavour non-singlet parts. The functions L2,fc describe the reactions where the 
virtual photon couples to the light quarks (u, d, s, n, d and s), whereas the if2,fc 
describe the reactions where the virtual photon couples to the c (c) quark. Hence 
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L2,fc and if2,fe in Eq. (1) are multiplied by e| and respectively. Moreover, when 
the reaction where the photon couples to the charm quark contains a light quark 
in the initial state, then it can only proceed via the exchange of a gluon in the t- 
channel. Therefore if2,g is a flavour singlet. This is in contrast with L2,g which has 
both flavour singlet and non-singlet contributions. Finally we have to add to Eq. (1) 
that part of the proton structure function F2{x, Q^) which contains charm quark loop 
contributions to the gluon self energies. The latter are only inserted in the matrix 
elements containing light partons (u,d,s,g). Combining this part of F2{x, Q^) with the 
contributions from L2^k leads to the correct asymptotic behaviour of the charm quark 
coefficient functions at ^ which allows us to perform the mass factorization 
discussed below. The NLO contributions to the charm quark coefficient functions 
were originally calculated in These functions were available analytically in the 
form of two-dimensional integrals which were then computed numerically. To speed 
up the integrations in Eq. (1) a two-dimensional grid of values for L2,fe and if2,fc 
together with an interpolation routine was provided in |Q. 

As a check on our exact formulae we have also calculated the asymptotic expres- 
sions for L2,fc and H2,k in the limit ^ m? 0. These asymptotic charm quark 
coefficient functions contain the large logarithms of the type \n^{Q'^ /rn^)\a' [ii^ /rn?) 
which dominate charm production far away from threshold. This has been shown in 
[^] by numerically comparing the asymptotic structure function defined by 

F™(3,x,Q^m2)= hm 

with the exact structure function in Eq. (1). Notice that in the last equation the 
exact charm quark coefficient functions are replaced by their asymptotic analogues 
mentioned above and Zmax = 1- The NLO results reveal that for Q"^ > 20 (GeV/c)^ 
and X < 0.01, i^^sYMP coincides with F2^^'^^ which implies that the large logarithms 
mentioned above entirely determine the charm component of the structure function. 
Since these corrections vitiate the perturbation series when gets large they should 
be resummed in all orders of perturbation theory. This procedure has been carried 
out in [0 and it consists of four steps. First we add the light parton (u,d,s,g) com- 
ponent of the proton structure function defined by F2{3,x,Q^) to i^A.SYMP -^^ 
(2). Second we apply mass factorization to the asymptotic charm quark coefficient 
functions containing the logarithmic terms ln*(Q^/m^) ln-'(/i^/m^). In this way the 
m^-dependence is separated from the Q^-dependence. The logarithms in the charm 
quark mass are put in the charm quark operator matrix elements whereas the log- 
arithms in the variable are transferred to the light parton coefficient functions 
denoted by C2,k {k = q,g). Both quantities have been calculated up to order in 
[|], and p respectively. In the third step we define new parton densities in a four 
flavour number scheme (FFNS) which can be written as convolutions of the original 
parton densities in the three flavour number scheme (TENS) with the charm quark 
operator matrix elements mentioned above. Therefore we obtain a new parton den- 
sity which represents the charm quark and is denoted by fc+c{z, fi"^). The latter has 
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the property that it does not vanish at /i = m in the MS-scheme contrary to what is 
usually assumed in the literature. In the fourth step we rearrange terms and obtain 
^ASYMP^3^^^g2^^2) ^ ^3(3, a;, Q^) = F2(4,a;,Q2) which is the FFNS result for the 

proton structure function. From the latter quantity one can extract the expression 
for the charm quark component of the proton structure function in the FFNS which 
will be denoted by 




where we have defined 

ClM = C^^linj) + nfClliuf) , Cl = rijCl . (4) 

The superscript PDF in Eq. (3) stands for parton density function which means 
that the charm component of the structure function can be completely expressed 
into parton densities multiplied by the light parton coefficient functions. Notice that 
F2^^ is a renormalization group invariant which implies that this structure function, 
like the one in Eq. (1) is explicitly independent of /i so that it becomes a physical 
quantity. Further it originates from the charm quark coefficient functions if2,fc since 
both are proportional to e^. The functions ^2,^, which are proportional to the light 
charge squared e^, also contribute to ^2(4, where they increase the number 

of flavours by one, and provide us with order matching relations between the 
TENS and FFNS light parton densities. The FFNS charm quark density is mainly 
determined by the size of the TENS gluon density G{3, z, fi"^). Therefore the latter 
plays a major role in the behaviour of F2^^^^ in Eq. (1) as well as -^2^^^ ^'i- 
(3). An analysis of both structure functions in 0] revealed that the former gives 
the best description of charm production in the threshold region where is small 
and X is large. On the other hand when is large and x is small it turns out 
that it is better to use F2^^ because it is in this region where the large logarithms 
ln*(Q^/m^) ln-^(/x^/m^) dominate so that they have to be resummed. Therefore the 
TENS is the most suitable scheme for the charm component of the structure function 
near threshold whereas far away from this region it turns out that the FFNS is more 
appropriate. 

One also needs a scheme which merges the advantages of these two pictures and 
provides us with good description of F2 c in the intermediate regime in Q^. This is 
given by the so called variable flavour number scheme (VENS). In [0 we proposed 
the following VENS structure function 

Q^ m') = Fl^^{4, x, Q') + F£^^^^(3, x, Q^ m') 

-F™(3,x,Q^m^). (5) 
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The above expression is a generahzation of Eq. (9) in 0, which was only presented 
in LO and has been implemented in a recent global parton density analysis |1T0| . 
(A different VFNS scheme has recently been proposed in 1T1[].) In and in |12| 
it was shown that F^™^ in LO is less sensitive to variations in the scale fi than 
each term on the right- hand- side of Eq. (5) separately. Further the VFNS scheme 
in LO has the properties that for ^ m^, F2^^'"^ — > _F^symp -^j^jc}! means that 

^VFNS ^ ^PDF ^i^iig q2 Q2 < ^2) ),ASYMP ^ pPDF Q ^^^^ pVFNS ^ 

F2^^^^- However the last relation is no longer true in higher order in The 
first reason is that all parton densities and the running coupling constant in -Ff^^ 
have four flavours whereas those appearing in ^2^^"^*"^ and ^MSymp j-^g^yg only three 
flavours. Notice that in lowest order we are allowed to put 0:^(3, yU^) = 0^(4, /i^) and 
G{3,z,fi'^) = G{4, z, fi"^). Hence in LO we can use the same number of flavours in 
all structure functions (PDF, EXACT and ASYMP). In NLO this is no longer true, 
which implies that there are regions in (small) where i^^sYMP ^gj-^not exactly 
cancel F^^^. The second reason is that new production mechanisms appear in NLO 
leading to the coefficient functions L2,k- The latter contain the prefactor el (e.g. 
the Compton process) and show up in F^^^*^"^ and _P^symp ^^^^ -Ff?^, which 
is proportional to el only (see Eq. (3)). Fortunately it turns out that both these 
corrections are small for < m? so that p^™^'^"^^ ~ ^exact,(2) j^^^^ierically. Yet 
another problem arises when one chooses a parton density set in which the charm 
quark density vanishes at /i = m. As we have found in |^ this property does not hold 
anymore beyond order in the MS-scheme. This can be remedied by either making 
a new parton density set satisfying the conditions presented in Eqs. (2.37)-(2.41) in 
[0] or by choosing a different scheme. As we use charm quark densities satisfying the 
condition /c+c(4, z^m?) = then we have to make an oversubtraction, which implies 
that the light flavour coefficient functions Cf ^ and Cf'^ in Eq. (3) differ from the 
results obtained in the MS-scheme. 

As an application we have studied the charm component of the proton structure 
function in the three schemes mentioned above. The coefficient functions used for 
these structure functions are all computed up to order (see |Q, 0, 0) so that 
we will denote them by ^^^'^^'^'(2) jT;q^ (^x), -^2^^^'''^^ in Eq. (3) and p^™^'^"^^ in 
Eq. (5). In Figs. 1-4 we have made plots of the x dependence of these functions 
for = 3, 12,45 and 170 (GeV/c)^ respectively and compared the results with the 
recent data from the HI p| and ZEUS [Q collaborations. We have chosen the next-to- 
leading log FFNS parametrization in |TH] with the condition that fc+c{z, /i^) = for 
/i < m for the input parton density set. Further we have taken the two-loop corrected 
running coupling constant with A4 = 200 MeV. Notice that in principle one has to 
choose a TENS parametrization for both the parton densities (see e.g Jl^) and the 
running coupling constant for the computation of F2^^'^^ and F:^^^^^. However we 
have checked that the latter are not significantly altered when we replace the parton 
densities in |]14| by those in |jl3|. Finally we adopt the mass factorization scale from 

■^This is only true if we put Zmax = 1 for F^^^^^ , which was not done in 0]. 
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P and use 

= + kQ\l-m'/Q^Y for > m\ 

= rn? for < rn? , (6) 

with k = 0.5, n = 2 and m = 1.5 (GeV/c^). 

At = 3 (GeV/c)^ (Fig. 1) the structure function -F™^'*'^'' is negative over the 
whole x-region so that we have not shown it. Nevertheless it is almost exactly equal 
to gQ i^j^ggg terms cancel in Eq. (5). Therefore at low the last two 

structure functions have no physical meaning because the parts in the heavy quark 
coefficient functions leading to the large logarithms do not dominate the QCD correc- 
tions to the structure function in Eq. (1). Hence at low the correct description is 
given by ^^^"^'^'^■(2) (^XFNS) which is almost equal to p^™^'^'^^ so that it is not nec- 
essary to use the VFNS description. However for larger Q^, i.e. > 12 (GeV/c)^, 
both ^^^^^"'"'(2) g^Yid ^2^°^'*-^^ are in agreement with most of the data. The same holds 

VFNS (2,') 

for which lies between the previous two structure functions. The only ex- 

ception is the value = 170 (GeV/c^) (see Fig. 4), where all descriptions are poor. 
The plots show the general features that p^^^'^^'^'^^ (TENS) is always below -^2^°^'*'^^ 
(EENS) which is due to the resummation of all the large logarithms in the latter 
structure function. Notice that these large terms are only included up to finite order 

. i;iEXACT,(2) 1 i;,ASYMP,(2) -r^. ,, , , r^VFNS,{2) r^EXACT,(2) , , 

m J^2,c ^iicl J^2,c ■ Fmally as expected ^ ~ -r2,c large x 

(threshold region) while p^™^'^"^^ ~ F^^^'^^-* at small x. However it is clear that one 
needs more precise data in finer bins of x and large to discriminate between the 
TENS and the EENS in order to observe the resummation effects incorporated into 
^2,^^'^"^^- The same holds for medium Q'^ where one would like to distinguish between 
the VENS and the other two schemes. 
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Figure Captions 

Fig. 1. F2™^'^'^\x, Q^) solid line and F2^^'^'^'^'^\x, Q^) dotted line plotted as func- 
tions of X at = 3 (GeV/c)^. The experimental point is from [Q. 

Fig. 2. F2™'(')(x, Q2) solid line, fI^^'^^\x, Q^) dashed line and F2™^'^')(x, Q^) 
dotted line plotted as functions of x at Q"^ = 12 (GeV/c)^. The experimental 
points are from [0[ closed circles and [^ open circles. 

Fig. 3. Same as Fig. 2 for = 45 {GeV/cf. 

Fig. 4. Same as Fig. 2 for = 170 (GeV/c)^. The experimental point is from 
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